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Vacuum Fluctuations of a Scalar Field during Inflation: 
Quantum versus Stochastic Analysis 

V. K. Oneml 

* Department of Physics, Istanbul Teehnical University, Maslak, Istanbul 34469, Turkey 

We consider an infrared truncated massless minimally coupled scalar field with a 
quartic self-interaction in the locally de Sitter background of an inflating universe. 
We compute the two-point correlation function of the scalar at one and two-loop 
order applying quantum field theory. The tree-order correlator at a fixed comoving 
separation (that is at increasing physical distance) freezes in to a nonzero value. 

At a fixed physical distance, it grows linearly with comoving time. The one-loop 
correlator, which is the dominant quantum correction, implies a negative temporal 
growth in the correlation function, at this order, at a fixed comoving separation and 
at a fixed physical distance. We also obtain quantitative results for variance in space 
and time of one and two-loop correlators and infer that the contrast between the 
vacuum expectation value and the variance becomes less pronounced when the loop 
corrections are included. Finally, we repeat the analysis of the model applying a 
stochastic field theory and reach the same conclusions. 

PACS numbers: 98.80.Cq, 04.62.-l-v 


I. INTRODUCTION 


MMC) scalar with a 
6|] several intriguing 


Quantum field theory (QFT) of the massless, minimally coupler 
quartic self-interaction on a locally de Sitter background predicts ll 
enhanced quantum effects of cosmological interest. A locally de Sitter background provides 
the simplest framework for an inflating spacetime. The cosmological constant A drives 
inflation in this framework and the scalar is just a spectator held. The operation of inhation 
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on the fluctuating quantum vacuum is to rip virtual particle pairs out of the vacuum by 
the Hubble flow before they find time to annihilate each other. After they emerge from 
the vacuum, they become real. As in flat spacetime, particles with smaller masses persist 
longer. Any massless particle which happens to emerge from the vacuum with sufficiently 
small wave number k<Ha{t), where H is the expansion rate, a{t) is the cosmic scale factor, 
and t is the comoving time, can persist forever 7-^. However, the particle production rate 
is suppressed for conformally invariant particles. Hence, this quantum effect during inflation 
is enhanced for the particles that are classically conformally non-invariant and effectively 
massless. The MMC scalars that we consider in this paper are an example of such particles. 


We computed 


the renormalized vacuum expectation value (VEV) of the stress- 
energy tensor for the MMC scalars endowed with a quartic self-interaction during inflation, 
at one and two-loop order. We showed that the renormalized energy density pren and pressure 
Pren of the scalar violate the classical weak energy condition p-|-p > 0 on cosmological scales 
at two-loop order. The equation of state parameter w = Pren/Pren < —1- As a result, a phase 
of super-acceleration is induced. This effect, however, is temporary since the field develops j4| 
a growing self-mass squared due to quantum processes and the particle production is cut off. 
The induced mass remains perturbatively small and does not go tachyonic which ensures 
the stability of the model . We studied how such a spectator scalar field might affect 
the measured curvature power spectrum A‘^{k,t) defined in terms of the Fourier transform 
of the equal time two-point correlation function of the derivative stripped three-curvature 
scalar field 7l(t, x) 

k^ 




More recently, we computed jo] the one and two-loop corrected power spectrum 

k^ 


27r2 




27r^ 


d'^xe * 0)|ff) , 


( 1 ) 


( 2 ) 


of the spectator scalar in our model. Correlations in various scalar potential models during 
inflation have been studied [n|. In this paper, we compute the one and two-loop corrected 
two-point correlation function of the spectator scalar 


{Q\ip{t, x)ip{t\ x')\Q) , 


(3) 


for t'<t and x' 7 ^ T, in our model. Bar on top of the scalar indicates that the ultraviolet 
modes with comoving wave number k>Ha are cut off the field. Hence, the IR modes with 
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H < k < Ha(t) are retained. The IR truncated field is guaranteed to reproduce the leading 
IR logarithms of scalar potential models in QFT. The advantage of using (p{t,x) is that it 
is a simple way of recovering the most important contributions; the disadvantage is that it 
leaves out sub-leading contributions. 

Physical wavelengths Aphys of the IR modes range from 27r/if to 27ia{t)/H. Had the 
lower limit of the k for the retained modes taken to be zero rather than H, the model would 
have suffered from IR divergences. To regulate these divergences we take the spacetime 
topology as T^~^xR where the coordinate toroidal radii R = 2tt/H. Here, the spatial index 
i = 1, 2,..., (D — 1). The VEV of the field (p(t,x) is zero at any event, hence the VEV of 
variation A(p{t, t'] x, x') = x) x') vanishes. We calculate the non-vanishing variance 


at one and two-loop order, using the two-point correlation function we compute in this paper. 
We also examine the stochastic properties of the IR truncated scalar in the same model 
repeating the computations in the context of stochastic field theory (SET) which is known 
to be successful in recovering the most important secular effects 1^ in inflationary QFT. 
Stochastic formulations of inflation have also proven useful in studying initial conditions 
|i^ . global structure 141, non-Gaussianitv 15|, power spectrum of a spectator scalar Q 
and conditions for eternal inflation 12], Starobinsky developed a stochastic formalism 
which reproduces the leading IR logs at each order in perturbation theory in scalar potential 
models. Starobinsky’s formalisms has recently been extended 19l-l21| to include scalars which 


interact with fermions 


22l | and with gauge particles 23|] but the inclusion of derivative 


interactions such as those of quantum gravity 24] and the nonlinear sigma model 2^, [26 1 
is a major unsolved problem. In this paper, we use the SFT approach of Ref. 27(| where 
the annihilation and creation operators in a quantum field are regarded as two complex 
conjugate random variables whose real and imaginary parts follow the Gaussian probability 
distribution with mean zero and standard deviation one. Results we achieved by applying 
the SFT are in perfect agreement with the ones we obtained applying the QFT. See also 
Refs 28|, l29(| cornpa ring different aspects of quantum and stochastic formalisms. 

A technique 30H33| to compute expectation values of functionals of a stochastic field 
is to integrate the functional weighted with a probability density function which obeys a 
Fokker-Planck equation. Employing this technique, we finally note a stochastic check of our 






















4 


two-point correlation function in the equal spacetime limit. 

This introduction is the hrst of hve sections. The outline of the remainder is as follows. 
In Sec. P we specify the background spacetime geometry and present the Lagrangian of 
the model. In Sec. uni we analyze the model applying the QFT: We compute the two-point 
correlation function and the variance of the IR truncated scalar held with a quartic self¬ 
interaction at one and two-loop order during inhation. In Sec. IIVI we analyze the model 
applying the SFT. We summarize our conclusions in Sec. |Vl The Appendices comprise the 
details of some computations in the paper. 


II. THE MODEL 


We consider a MMC scalar with a self-interaction potential V(ip). The Lagrangian density 
of the model is 

£ = -^(l + dZ)d^ipd^ipg'^’'^/^-V((p)x/^ , (5) 

where 6Z is the held strength counterterm and g/j,iy is the background metric of locally 
de Sitter spacetime. The invariant line element can be expressed in comoving coordinates as 

ds^ = gfj_iydx^dx‘' = —dt^ + a^(t)dx ■ dx , (6) 


where the scale factor 

a(t) = , 


(7) 


is normalized to unity at the initial time ti = 0. We work in a D-dimensional spacetime, 
hence the indices fi, 1 ^ = 0,1,2,, (D — 1). Thus, x^ = (x^, x), x^ = t, and = (9o, V) in our 
notation. We use the open conformal coordinate patch of the full de Sitter manifold where 
each of the (D — 1) spatial comoving coordinates x* lies in the range 0<x^ <H~^. 

The equation of motion for the scalar held 


V3(t, x) + (D-l)Hip(t, x) -^ip(t, f) + ^ ^ ^ ’ 

1+dZ 

where a dot denotes derivative with respect to comoving time t, has the solution 


ip(t,x) = ipo(t,x)— I dt'a^ ^(t')jd^ ^x'G(t, x;t', x') 


,, ^,.V’(p^)(t’,x’) 


1 + 6Z 


( 8 ) 


(9) 


The free held ipo(t,x) in Eq. ([9]) obeys the linearized equations of motion and agrees with 
the full held at f = L = 0. The Green’s function G{t,x-,t',x'), on the other hand, is any 
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solution of the equation 

V2 

Git, x; t', x') + (D — l)HG(t, x; t', x') -^0 =5{t—t')5^~^{x—x') , (10) 

which obeys the retarded boundary conditions. We continue studying the model applying 
the QFT in the next section. We study it applying the SFT in Sec. [IVl 

III. QUANTUM FIELD THEORETICAL ANALYSIS 

In this section, we analyze the model described in Sec. |TT] in the QFT point of view. 
Section IIII Al is reserved for the non-interacting limit of the model. Section IIIIBI considers 
the model with a quartic self-interaction potential. 


A. Free Theory 


Recall that the background is homogeneous and isotropic with a topology ^xR and 
we work on a comoving coordinate patch where the D — 1 spatial coordinates x^ he in the 


hnite range 0 <x^ < H The modes of helds are, therefore, discrete as they are for 
hnite spatial manifold. Hence, we expand the free scalar held in a spatial Fourier series 

P 




d. 


(po{t,x)=H 2 W- 


{D-l)Ha^-^{t) 


n^O 


+5^ Ut, +u*it, \ , (11) 


where the comoving wave vector k = 2TiHn and the comoving wave number k = ||/c||. The 
mode with physical wave number /cphys(t) = k/a{t) crosses the horizon for the hrst time at 
t = tk- Hence, its physical wavelength Xphys{tk) = ‘2'Kk~^y^{tk) = H~^{tk)- The expansion rate 
H, however, is constant during de Sitter inhation, H{tk) = H. Therefore, this mode has 
k = 2TTHa{tk) = 2TTHe^^'‘. The Bunch-Davies mode function u{t,k) in Eq. flTTl) is a solution 
of the spatial Fourier transform of the linearized held equation. 


^o{t, k) + (p-l)H^o{t, k)- 




a2(t) 


(po{t,k) = 0 , 


( 12 ) 


for the modes with k^O. It can be given in terms of the Hankel function of the hrst kind as 

k 


u{t, k) = i 


TT 


AHa^-^t) 


n 


( 1 ) 

D-l 

2 


Ha(t) 


(13) 
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For the fc = 0 mode, the two solutions of Eq. flT^ are a constant and Hence, the hrst 

two terms in Eq. (1TT|1 represent the zero mode of the quantum held (po(t, F). The operators 
Q and P satisfy the commutation relation 

lQ,P]^i- (14) 


They can be expressed in terms of the Fourier transforms of the held and its hrst time deriva¬ 
tive for k = 0, evaluated at f = 0, 

0o(O,O) + p-l)i7^o(O,O) 


Q = 




P = 


(Po(0,0) 


(15) 

(16) 


The annihilation and creation operators, and At in Eq. fllip . may change as diherent 
Green’s functions are used in Eq. ([9]), but their nonzero commutation relation 


-4a. -4t 




(17) 


remains hxed. They can also be expressed in terms of ipo{0,k) and (pQ{0,k). Recall that 
Eq. flTT)) and the canonical commutator 


(po{t,k),(pQ{t,k') 

hxes the Wronskian 

u{t, k)u*{t, k) k)u{t, k) = 

Using Eqs. flTT]) and flT^ we obtain 




7.^-1 m 


A- — 

—t 



A^ — 

i 


77 V ^ 


'U*(0, /c)</9o(0, k) —M*(0, /c)0o(O, k) 
^(O, fc)^o(0, —fc)—n(0, /i:)(po(0, —k) 


1 



(18) 


(19) 


( 20 ) 


The state |G) annihilated by all Aft is known as Bunch-Davies vacuum. What we assume 
about the dependence of vacuum state |f2) on the zero mode operators Q and P does not 
matter much because there is only one zero-mode as opposed to an ever-increasing number 
of nonzero modes. Thus, we neglect them in free held expansion flTTl) . 

Free held expansion flTT]) also includes arbitrarily large wave numbers which cause ultra¬ 
violet divergences. Therefore, we may consider the expansion obtained by cutting out the 
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sub-horizon (ultraviolet) modes and retaining only the super-horizon (infrared) ones 


n^O 


( 21 ) 


In the above expansion, we may further take the IR limit of the mode function flT^ given 
in terms of the Hankel function 


«)t.( 2 ) = ^2=1(2)+»Wz) 


( 22 ) 


Here z = jArpr, J p-i is the Bessel function of order and 


D-l 


Ha{t) ■ — 


Yd-Az) = Jp-iiz) cot ( --- ) — J (p-i\(z) CSC 


(23) 


is the Neumann function of order Using the series expansion of the Bessel function 


Jp-i(z) = 


z\ — 


D-l 


^1) 


+ 0{z^ 


the reflection formula for the Gamma function 


(24) 


pvTcsc(p7r) =r(l—p)r(l-l-p) , 


(25) 


with p=^-^ and the identity 


D-l\ ^ r(D-i) 


2 ) 2^-2 r(f) 


D\ 1 


one gets 


22 


23 I. l3l| the leading IR limit of the mode function 


M(t, k) 


r(D-i) H- 


l + o 




(26) 


(27) 


r(f) ( 2 /t)l°->)/p‘'''Vfl- 2 a 2 (() 7 | ' 

Then, substituting u{t, k) in Eq. flTT)) into Eq. fl^ we obtain the IR truncated free held in 
D-dimensions, 


■pi'T^ 1 5 tjD—3/2 , 




(28) 


Let us dehne the hrst time derivative of the held as 




r(D-i) H 


D-h 


rfR') 2(-^“b/2/t(i^-3)/2 




(29) 





























where and n = ||n||. Its stochastic analog represents stochastic noise in Sec. II V Cl 

We now want to calcnlate the two-point correlation fnnction of the IR trnncated free field 


{n\Lpo{t,x)^po{t',x')\n), 


(30) 


for t'<t and x'^x. It will be needed in Sec. IIIIBl where we compnte the correlation fnnction 
of the IR trnncated fnll field. Using Eqs. (ITTll and (EHD, we obtain 


{Q\ipo{t, x)ipo{t', x')\Q) 


m^O 


r2(D) 20-1 ^ A,(D-l)/2 


5^' 

n^O 


y / -ik'-x'x /Qi \ 


Here k' = 2nHm. When m = nwe have k' = k, therefore, 


{n\ipo{t,x)(po{t\x')\n) 

r2(D-l) if ^-2 


r2(f) 220-%^- 


tE 


Q{Ha{t)-H27in) Q(Ha{t')-H27m) j 2 M(x-x') 


n^O 


n 


D-1 


. (32) 


The discrete snm over n in Eq, fl52]) can be approximated as a (D —l)-dimensional integral 


E 


eiHait)-H27rn) Q{Ha{t')-H27:n) ,,^Hn.(x-x') 

nD-i 


^dk 


— 0(iia(t)-A;)0(iia(t')-fc)U^^^“®(^) , (33) 


k 


where Ax =|| Ax || = || x —x'|| and 9 is the angle between the vectors n and Ax. Evalnating 
the angnlar integrations on the right side of Eq. (13^ yields 


D 1 ^1 r(f) 


'dk^,^^ ,,sin(/cAx) 

—Q{Ha{t) — k) Q{Ha{t ) — k)- 


r(D — 1 ) 7 q k '' '' ' ’ '' '' ' ' kAx 

Becanse 0<t'<f, we have Ha{t')<Ha{t) and hence, 

0 {Ha{t)-k)e{Ha{t’)-k) = e{Ha{t’)-k) . 


(34) 


(35) 


Therefore, the integral over k in Eq. 


dk . sin(fAx) f^°‘^^'^dksm{kAx) 

0 (iia(t ) — k)- — ' — ' 


k 


kAx 


' H 


k kAx 


./ N sin (a(t')h/Ax) , 

= c.(a(OJ?A.)--c.(J?A.) 


HAx y y 

sin(iiAx) 


HAx 


(36) 




















9 


where the cosine integral 

^cos(t) 


/ “ cnsiti —1 ^ ^ f — 

dt -^-= 7+ln(Q)+^ 2n(2n)! ’ 


t 


is an entire function. The Euler’s constant 7 ~ .577. Employing Eqs. 
Eq. fl3Tll we obtain 


and 


E 

n^O 


e{Ha{t)-H 27 m)Q{Ha{t')-H 2 Tin) ^ 


i2'KHn-{x—x') 


D-l 


n 


r(f) 


~ 2"“‘7rT 


rp-i)L 


./ X sm(a(t')HAx) , smiHAx) 

ci{a{t’)HAx) -- --ciiHAx)+ ^ ^ 


a{t')HAx 


HAx 


We can power expand the time dependent terms in Eq. 

n=l ^ ^ 

with a = a{t')HAx. Substituting Eq. (j39ll into Eq. (]38|1 yields 

E 

n^O 


e{Ha{t)-H 27 rn) e{Ha{t')-H 27 rn) ^^Hn.(x-x') 


n 


D-l 




r(D-l) 


C(Aa;)+ln(a(T))+y^- 


2n 


n=l 


2n ( 2 n+l)! 


where we dehne the comoving spatial separation dependent function 

siniHAx) 


C{Ax) ='Y — l+\n{HAx) — ci{HAx) + 


HAx 


■E 




i2n 


2 n( 2 n+l)! 


n=l 


(37) 


m 


(38) 


(39) 


(40) 


(41) 

(42) 


Combining Eq. (13^ and the analytic form of the sum given in Eq. (I38ll gives the two-point 
correlation function of the IR truncated free held as 


HD-2 T{D-1) 


2 ^-i 7 rW 2 r(f) 


, sin (a(f')-f^Aa:) , sm(HAx) 

ci(a(t’)HAx) - ; \ \ - ^--cUHAx) + —A —L 

TTA HAx 


(43) 


a{t')HAx 

We can use this analytic form of the correlator to infer its behavior for a fixed comoving 
separation (increasing physical distance) at late t' when a{t') becomes large. Equation 
and the asymptotic form of the function ci(a) for large argument 

sinfa) 


sin(a)/ 2! 4! 

ci(a) =- 1 - - + — 

a \ 0(2 0(4 


cos(q() f 1 3! ^ 5! 


a \aa^a^ 


a 


(44) 




























10 


imply that the correlator at a fixed comoving separation freezes in to a nonzero constant 
value at late times. The usual result in QFT is for it to go to zero as the physical separation 
increases. 

Combining Eq. (HOjl and Eq. (H2ll in Eq. (13^ gives the correlation function as a series 
expansion 




H 


D-2 


T{D-1) 


2D-1^D/2 

hd -2 r(D-i) 


C(Aa:) + ln(a(T))+y^- 


n=l 


2n (2n+l)! 


2D-1t,DI2 




(45) 
. (46) 


2 n(2n+l)! 

^=l ^ ' 

Equation fl^ for Eq. fHOjl l can be used to interpret the correlator at a. fixed physical distance. 
As in Ref. (27|, we may choose Ax so as to keep the physical distance a{t')Ax a constant 
fraction K of the Hubble length, 


Ax = 


K 


Ha{t') ’ 


(47) 


and use Eq. fl39l) in Eq. fjTHl) to obtain 

{n\(po{fix)ipo{t\x’)\n) ln(a(F))-ln(A:) + ci(A:)- ^^^f^^ -7 + 1 


K 


-E 




r(D-i) 


2D-i^d/2 


ln(a(t'))+E 


n=l 

^_l)n^ 2 n 

'-^2n (2n+l)! 


2n (2n+l)! 


(l-a-“( 0 ) 


(48) 


(49) 


Notice that the time dependent part in the infinite sum of Eq. fITHD (or Eq. fl49|) ) decays 
rapidly as 0{a~‘^{t')). Since ln(a(f')) = Ht', we infer that the correlator at fixed physical 
distance grows linearly with comoving time. This is because more and more particles are 
created, which increases the local field strength. At this point, let us note that the equal 
space limit of Eq. fHbll yields 


{n\ipo{t,x)(poifi f) IH) = JimJH \(po{t, x)ipo{fix')\9.) 

X '—>-X 


hd -2 r(D-i) 
-2D-It,d/2 y(A) 


ln(a(f')) , 


(50) 


and taking the equal time limit of Eq. fl50|l gives 


T{D-1) 


(H|y?o(^)^)l^)=jim(f^|(^o(C^)<^o(^',^)|^) - 2 D-i^d/2 • (51) 


The stochastic analog of the two-point correlation function for the fcpfit, x) will be needed 
in Sec. IIV Cl For the sake of completeness, let us compute it here, in the context of QFT, 
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as well. Equations ffT7|) and fl29|) imply that 

r2(E)-l)i72D-i 




n^Om^O 


^Ha(tn),k ^Ha(t'^),k' ik-x „—ik'-x' 


p2(B) 2^-1 f^'iD-3)/2 


■e e 


^n,rn j (52) 


where t'<t, x'^x and k' = 27THm. When rn = n we have k' = k, therefore, 




_ ^ ^ ^^Ha(tn), Ha{t'^) ^i27rHn-{x—x^) 


n 


D-3 


p2(D) 220-4vr^-3. 

We make the integral approximation to the discrete mode sum in Eq. (15^ to obtain 


(53) 


E 

nj^O 


^ H a{tn) , H ^i2nHn-{x—x') 


n 


D-3 




^ dkk 
{2'KHy 


5{Ha{t) — k) 5{Ha{t ) — k)e 


— cos ( 61 ) 


(54) 


Evaluating the angular integrations on the right side of Eq. fl5^ yields 


2d^d/2 r(|) 


dk 6(Ha{t) — k) 6{Ha{t') — k) 


sin(/cAa;) 

Ax 


(27r)3i72r(Zl-l) 

Then we perform the remaining integral over k in Eq. (]55ll . 

/ dk6{Ha{t) — k)6{Ha{t) — k) --- = --- -5{Ha{t) —Ha{t)) 


(56) 


Ax 


where we used 5{Ha{t) — Ha{t')) = Using Eqs. fl5^ and fl56|) in Eq. 

End 


Ax 

sm{a{t')HAx) — 
a{t')HAx ir~ ’ 


(56) 


, we 


i 2 M(x-x') _ 2^7r^/^ r(f) sin(a(tOijAa:) 

^ n^-3 {2ttH)^T{D- 1) a{t')HAx ^ 


(57) 


Finally, inserting Eq. flFT)) into Eq. fl5^ gives the correlation function we wanted to compute. 


The equal space limit of Eq. fl58|) yields 

{^\Uo{t^x)f^^{t',x)\Q)c:^^^j^^Y^^6{t-t') . (59) 

This brings our discussion on the free theory to an end. In the next section, we study the 
interacting theory with a quartic self-interaction potential. 
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B. Interacting Theory 


The solution of the equation of motion for a MMC scalar with a self-interaction potential 
V{ip) is given in Eq. ([9]) as 

(p{t,x) = ipoit,x)-J d^~'^x’G{t, x; t\x') ^ ^ , 


where the Green’s function, defined in Eq. (fTOj) . is 19j 


G(t, x; t\ x') = — (po{t, T), (po{t\ x') 


( 60 ) 


To get the IR truncated full field ip{t,x), we need the leading order IR limit of the 
G{t,x]f,x'). Combining Eqs. flT^ . fl2^ and fl23ll we hnd 

TT I 

u{t, k)u*{t\ k)—u*{t, k)u{t\ k) =i- 


71 CSC^ 2 


H ]D-1 


2H it) it') 

1 / k^ 


J D-l (z) J — (z 


k^ 


ma^ity H‘^a^it'] 




, (61) 
• (62) 


Equations (l2Tll and (|25ll have been used in obtaining the limit in the last line. Hence, 


plugging Eq. fl6^ into Eq. fl60|) . retaining only the leading terms, gives 


22 


23 


3l| 


G{t, x; T, x' 


— ^{x—x') 


iD-l)H 




(63) 


When this limit is inserted into Eq. (E]) it is multiplied by the integration measure 
Then, the hrst term in the square brackets, = 1, contributes over the whole 

range of the integration. The second term proportional to on the other 

hand, contributes signihcantly only for t' ~ f, and hence, negligible in the leading logarithm 
approximation we consider. Thus, Eq. ([9]) yields the IR truncated full field as 

1 


<p{t,x) = ipo{t,x) 


_ v'm\s)) 

{D-1)hJo 1+dZ 


(64) 


We consider a scalar held with quartic self-interaction, hence 


( 65 ) 
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where A is the coupling constant, and are the coupling constant and the mass 
squared counterterms, respectively. Counterterms cannot contribute in leading logarithm 
order. To see this, let us hrst note that bZ ~ C>(A^), bvr? ~ C^(A) and hA ~ C>(A^) in the 
model jj]. In the ^(A+hA)^?"^ term, A and 5A multiply the same Therefore, contributions 
involving and contributions involving must have the same IR logarithm structure. 
Because the latter are suppressed by at least one extra factor of A, with no more IR log¬ 
arithm, they cannot be in leading order. Next, let us compare the and terms. 

Although, 5m^~A, it multiplies p‘^. Therefore, at a given order in A, contributions involv¬ 
ing the former cannot produce as high order IR logarithms as contributions involving the 
latter. Finally, the held strength counterterm 5Z appears in the held equations in the form 
V{p)/ {1 + 5Z) = V'{p)\l — 5Z + {5Z)‘^ — ■ ■ ■] where hZ~(9(A^). Consequently, exactly the 
same leading logarithm contributions is obtained from the simplihed held equation without 
the counterterms. 


p{t,x)=po{t,x)- 


A 


Q{D-l)Hj 


dt p'^(t'\x) . 


( 66 ) 


Iterating Eq. (l66ll twice, we hnd the IR truncated full held in terms of p^it^x) as 


p{t,x) = pQ{t,x) 


I 


6{D-1)HJ, 

'24(T)-l)3if3 


dt'pQ{t\x) + 




12(D-l)2ifVo 

rt' 


dt'pQ (F, x) / dt"pQ (t'\ X 


dt'p^o{t\x) dt"p^oit'\x) dt'"p^oit"\x) 


+ 2 / dt'po{t',x) 


dt''pl{t", x) 


+ 0(1^). (67) 


Using Eq. fl67D . we can compute the two-point correlation function of the IR truncated full 
held. Up to 0{l^) we have 


(ff |(p(f, x)p{t', T') l^d) = (f^|7^o(C x)po{t', f') 1^) 


A 


6{D-1)H 

A2 


(f2|(^o(U ^)[ dtpl{t,x')\Q,) + {Q.\ f dt''pl{t'\x)po{t\x')\Q,) 
Jo Jo 

(r2|(po(t,x) f dipl{i,x') [ dip\{i,x')\Vt) 


12(0-1)2772 


-|-(r2| / dt"pQ{t'\x) dt"'pQ{t'",x)po{t\x')\^l) + -{Q.\ dt"pQ{t",x) dtpQ{t,x')\^l) 


, ( 68 ) 


where 0 < t'" <t" <t, 0<t <t <t' and t' < t. The perturbation theory breaks down when 
ln(a(t)) =77f ~ 1/\/A jl|. So, for A -Cl, one can have a long period of time during which the 
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perturbation theory and hence Eq. (168|) are valid. In the next three sections, we compute 
the tree, one and two-loop contributions in correlation function fl68|) . 


1. Tree-order contribution 


The tree-order contribution in the two-point correlation function of the IR truncated full 
held is just the two-point correlation function of the IR truncated free held which is obtained 
in Eqs. (145]) and (145]) as 


HD-2 y{D-1) 

-2D-It,DI2 rg) 

hd -2 r(D-i) 

-2D-It,DI2 y(R) 


{Q\(po{t,x)ipo{t',x')\Q) 

./ / /X X Sm(a(t')HAx) , 

ci{a{t)HAx) -———- ci{HAx)T 

a[v)HAx 

n=l ^ ' 


sin(/fAa;) 

HAx 


(69) 

(70) 


Quantum correcting this classical result, however, requires some work. In the following two 
sections we carry out the calculations for the one and two-loop contributions, respectively. 


2. One-loop contribution 


Computing the one-loop contribution to the tree-order two-point correlation function in 
Eq. (168]) involves evaluations of two VEVs. The hrst one is 

{^l\(fo{t,x) f di(fl(t,x')\^l)= f di{n\(fo{t,x)(pl{i,x')\^l) . (71) 

Jo Jo 

It can be calculated either by directly inserting free held (l28]) for each power of <^o in Eq- (TH!]) 
and using the commutator algebra (ITT]) or, equivalently, by topologically inequivalent field- 
pairings which yield 

f di (fl|(,7o(C x) ^01^) — [ diS ■l{Q\(po{t, x) <^o(C ^0 

Jo Jo 

'-l)^{a{i)HAxy^ 




2‘iD-2^D 


r^(f) 


'0 


C (Ax)-h In (a(f))^ 


n=l 


2n (2n-f 1)! 


ln(a(f)) , (72) 
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where we used Eq. fHSjl . the fact that t <t and Eq. fl5T|) . Evaluating the integral in Eq. fl72|) 
we obtain 


dt x) x')\il) = 


H 


2D-5 f2 


r^(E)-l)^j21n^(a(t')) 


220-%^ r2(f) 


+C(Aa;)ln"(a(t')) 




[-l)^(HAx) 


2n 




2^ n2(2n + l)! 

The second VEV which contributes at one-loop order in Eq 
similarly, 




(73) 

can be calculated 


(r2| / dt"(pQ{t",x)ipo{t\x')\Q.)= / dt" 3 ■l{^l\ipo{t'\x)(po{t\x')\^l){^l\ipQ{t",x)\Q.) 

Jo Jo 

where we used Eq. fISTD in the last equality. We have 0 <t" <t and t' < t. Therefore, in 
order to evaluate the remaining VEV in the integrand we break up the integral into two 
as ^^dt” = Jq dt" + J^idt”. In the first integral t” < t', whereas in the second t' < t". Hence, 
appropriate use of Eq. (H5l) in Eq. gives 

{n\f dt"(fl{t'\x)(po{t',x')\n) 


H2D-4 r2(T)-l) 

22D-2^D ^2(0) 


3{C{Ax) dt"ln{a{t"))+ dt"ln\a{f)) + ln{a{t')) dfln{a{f)) 


'0 
\ 2n 


l^{-ir{HAxr 


2^ n{2n+l)\ 

71=1 ^ ^ 

We evaluate the integrals and hnd 


dtV’^(t")ln(a(t")) + a"’^(t') / dt"ln(a(t")) 


(75) 


7' /72D-5 p2m_i\ 

(H| / dff^{t'\ x)(po(t', X') |f^) = 


22D-1t^d 


C(A)r)-|-ln(a(T)) ln^(a(t))- 


ln^(a(t')) 


2^ n(2n+l)\ 

n=l ' ' 


i2-(t')fln2(a(t))-ln2(a(t')) + ^^^^^^ ^ 




2n2/ 2n2 


(76) 


Combining Eqs. fH2]) . (175]) and fl75]) in Eq. (165]) yields the following total C>(A)-correction 


{Q\(p{t, x)ip{J, X ') |f^)l-loop ^ ■ 


A \hA{a{t')) . 


(L)-l)22^7r^ r2(A) 


-|-ln(a(T))ln2(a(t)) 


2^ n(2n+l)\ 

71 = 1 ^ ' 


2n 




ln2(a(t))-ln"(a(T))4 


-hV(a(t))-ln2(a(t')) 


,,, 21n(a(T)) 1 


n 






(77) 
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As in Eq. (H3|) . to infer the result at a fixed comoving separation at late times during inflation 
we must fix Ax and then take t' large. (Recall t>t'.) The asymptotic form of the one-loop 
correlator is obtained as 

A y^{D- 1) f 41n^(a(f )) 


(f2j(p(t, x)ip(t', f') l^)l-loop ^ ■ 


p_l)22^7r^ r2(f) 


ci(HAx) 


sin (HAx) 
HAx 


ln^(a(f))-|-ln^(a(f')) +2 ln(ifAa ;)-|-7 —1 ln^(a(f')) 


a 


a‘ 


a 


— 3 J -4 1, 1 , 1; 2, 2, 2, ln(a(t')) + — 4 -T 5 1, 1 , 1 , 1; 2, 2, 2, 2, — 


a‘ 


12 

H'^Ax^ ( 5 H’^Ax^ 

1,1.1.1; 2.2,2.2.-; 


(78) 


where we used the large argument expansion of the ci(a) function with a = a{t')HAx. 
The mJ^nS are generalized hypergeometric functions. Equation (175]) implies that one loop 
correlator grows negatively for a hxed comoving separation during inflation. 

To obtain the one-loop correlator for a hxed physical distance, we choose a{t')Ax a 
constant fraction K of the Hubble length as in Eq. flT7|) . 

A 21n3(a(f')) 


sin (A) 


ci(A)- 


K 


ci{Ka 


(T)-l)22^7r^ r2(f) 

sin(Aa“^(t')) 


Aa-i(t') 


ln"(a(f)) 


-pi(A')-^^-21n(A-) + 


ci(Aa 


sin(Aa 


!-27l 


+ 2-27|ln2(a(t')) 


/ 5 KA A^ / 5 A^ 

37-4 R, 1 , 1 ; 2, 2, 2, j ln(a(t')) + — 4 -T 5 (l, 1,1,1; 2, 2, 2, 2, — 


6 


A2 / 5 

-4^5 1,1,1,1; 2, 2, 2, 2,-; 


2' 4 
A2 


12a2(t') V’ ’’’’’’ ’ 2 ’ Aa?{t') 

We employ Eq. (15^ in this analytic form of the one-loop correlator and get 

A A2^-6r2(D-l) J ln^(a(f')) 


(79) 


{VL\if(t, x)ip{t\ X') |H)i_ioop ^ • 


(D-l)22^7r^ r2(f) 


f ln(a(f'))ln^(a(f)) 


°° /_ 1 \ nT^2n °° /_.. \nT^2n °° ( _i \ nT^2n 

+ rin2(a(t))-ln2(a(t'))l V-M-^ + ln(a(T)) -rr - V-M-rr 

L V V V V 2n+l ! ^ ^ ’’^n^{2n+l)\ ^2n^{2n+l)\ 

n=l ^ ' n=l ^ ' n=l ^ ' 

where the terms in the last line decay exponentially during inhation. Equation fISOD tells us 
that the one-loop correlator at hxed physical distance grows negatively during inhation. 
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In the next section, we compute the two-loop contribution to the correlation function of 
the IR truncated scalar held. 


3. Two-loop contribution 


Computing the two-loop contribution in two-point correlation function 
nations of three VEVs, of which the hrst is 

{Q\(pQ{t,x) [di(pl{i,x') [ di(pl{i,x')\Q)= fdifdi< 

Jo Jo Jo Jo 

2-3-1 {n\(po{t,x)ipo{i,x')\n){n\Lpo{i,x')ipo{i,x')\n){n\Lpl{i,x')\n) 

-h3-2-l {n\(po{t, x)(po{i, x')\^l) {il\cpo{i,x')(po{i,x')\il) 


involves eval- 


( 81 ) 


where we have 0 < t < t < t' and t' < t. This time ordering makes it easy to read off the 
VEVs in the integrand from Eqs. (H5|) . (I50D and (ISTD . Substituting the outcomes yields 

rt' 


2n 


I ~ ~ ~ ~ pjoD-o r^fD—li 

(f2|^o(C x)dix')jdi^l{i, f') 1^) = 200 - 3 ^ 30/2 

IdtJdh-S-l C(Ax) + l„(a(f))+|;VElg^ 

pt' pi _ r 00 

di dt 3 ■ 1 ■ 1 
Jo Jo 

f di fdi3■2-1 
Jo Jo 


C(Aa;)-|-ln(a(t)) 


2n 


n=l 

00 


C(Aa;)-l-ln(a(f)) 


n=l 


-l)^{a{i)HAx) 
2?7,(2n-|-l)! 

2?7,(2n-|-l)! 


ln^(a(f)) 
ln(a(t))ln(a(t)) 
ln^(a(t)) 


(82) 


The results of the integrals in Eq. 
and (1A311 in Eq. (|8^ we obtain 


are 


given in Appendix I A II Using Eqs. flAip . 


(f2|(^o(C x)f di (pl{t,x') f di (fl{i, T') |f2) 

Jo Jo 

H30-8 y^{D- 1 ) [91n'^(a(t')) llC(Aa ;)^^^4 ^{-l)^{HAxy^ 


230-2^30/2 Y^{—) 




ln3(a(0)4 


31n^(a(T)) 211n(a(T)) 33 


4n 


8n? 


16n3 


^ n2(2n-|-l)! 

T.= l ^ ' 

^ 31n^(a(d)) 31n(a(d)) 33 


8n 


2n2 


16n3 


























18 


The second VEV which contributes at two-loop order in correlation function 
evaluated similarly, 


IS 


rt ft" 

-t\ / T.w-S/'j.W 


i ft" 

" / Jj-t" 


(fi| / dfiplit”, x) / dfiflit”', x)ifo{t\ f') 1^) = / df / dt 

Jo Jo Jo Jo 

2■ 3■ 1 (flI(po(t" x)ipo{t\x')\^l) {^l\ipo (t", x)(fo (t'",f)Ifl) (flI(po(t'",f)Ifl) 
-hl-3-1 T)|fl) {il\(po{t''\ x)ipo{t\ x')\il) f)|n) 


+3-2-1 


mifoit”, x)^po{t''', x)\il) {n\cpo{t''',x)cpo{t',x')\n)> , (84) 


where 0<t'"<t"<t and t'<t. Employing Eqs. (|5nil and (|5Tll in Eq. (|8T|l we get 

ft ft" tt 2D—A Y‘^(T) _ 1 

(f2| / dt"+o(^" / dt"'f^{t"\ f)+o(f', T') 1^) - 


'0 


'0 


22D-2^D 


-t ft" 


/ dt" dr2-3-l(f2|+o(t" x)+o(t',x')|f^)hi^(a(r)) 

Jo Jo 

Jo Jo 

~^Jdt"J dt"'?)-2-l{Vt\(po{t"',x)(po{t',x')\VL)\Y?{a{t"')'\ . 


(85) 


Because the ordering between t' and t" and between t' and t'" are unknown in the double 
integrals, we cannot just read off the VEVs from Eq. flTSD . We can, however, decompose the 
double integrals so that a dehnite ordering between the time parameters exists. That makes 
the evaluation of the integrals, using Eq. (IT5|) . possible; see Appendix lA 21 Substituting 
Eqs. (1A6I1 . (lAllI) and (1A13|1 into Eq. (l85il . we find 


/r^i /■*,// -2. // ,n-Z, m T^{D-1) [llln®(a(f')) 

{^l\ dt , x) dt ipl{t ,x)(po{t,x)\^l)^ 

Jo Jo 

,3/ 




— ln'^(a(f'))ln(a(f)) ■ 


-ln^(a(f))4 


23D-2^3D/2 r3(|) Y 

ll[ln(a(f'))+C(Ax)]. 


20 


2 ' 4 

1 / /\\ 1 4 / / / / NN 51n^(a(f')) llln"^(a(f)) 

21 n®(a(T)) -21n^(a(f'))ln(a(f))- 


'h r'' 

31n^(a(f')) , 2/ 1 Alln^(a(t')) ou.2, 31n(a(f')).2 


ln^(a(f)) — 
4 n 




-3hr(a(T))ln(a(f))- 
, 2 / 


ln^(a(f)) 


1 f33W{a{t')) 31n^(a(t))5 |_1 | ^331n(a(t0) 31n(a(t)) 5| ^ 33 


8 


8 J n^\ 8 2 

31n^(a(f)) 31n(a(f)) 33 


lOn^ 


8n‘^ 


2n3 IQn^ 


( 86 ) 
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Finally, the third VEV which contributes at two-loop order in correlation function fl68|) is 


pt pt/ pt pt/ 

{Q.\ dt"x) diipl{i,x')\^l)= / dt” di< 

Jo Jo Jo Jo 

3-1-1 {^l\ipo{t", x)(fo{i, x')\Q) (n|<^o(t",f)|n) f')l^) 

-h2-3-l (n|<^o(t",x)|fl) {Q\(po{t'\ x)(fo{i, x')\^l) (fl|<Po(^f')|n) 


3-2-1 


(f2|<y9o(t",x)<y9o(t,f')l^) 


(87) 


where 0<t"<t, 0<t<F, and t' <t. The integrand on the right of Eq. flHTI) consists of three 
terms that are added together. Notice that the second term is twice the hrst term. Using 
Eq. (l5Tll in Eq. (|87ll we get 

(121 [dt"ipl{t",x) f di(pl{i,x')\Q,) 

Jo Jo 

/ t |•t' f j^2D—4: pS/n 1 \ 

dt"j^ di I > ln{a{t"))ln{a{i)) 


+6 


{Q\ipo{t",x)ipo{t,x')\Q) 


1 3 


( 88 ) 


As in the evaluation of the second VEV given in Eq. fIST)) . we need an ordering of the time 
parameters t” and t to evaluate the expectation values in the integrand on the right of 
Eq. fl88|l . Therefore, we similarly break up the double integral into two parts: JgdF'/Q di = 
f^dt"f^ di+ J^dt"jQ di where in the hrst part on the right t" < f, whereas i< t" in the second. 
Hence, employing Eq. flTSl) appropriately, we obtain 

H^d-o Y^iD-1) 


(121 / dt"(pl{t",x) f dt(pl{t,x')\Q,) 

Jo Jo 


j dt"J dt{9 


C(Ax) + ln(a(2"))+X] 


23D-3^3D/2 r3(':P^ 

-l)^{a{t'')HAxf^ 

n=l 


+6 


C(Ax)-hln(a(2"))+^ 


2n (2?7,-|-1)! 

-l)^{a{t")HAxf^ 


ln(a(2")) ln(a(2)) 

3- 


n=l 


2n (2?7,-|-1)! 


dt" j dt{9 


C(Aa:)-Hln(a(t))-[-y^- 


-lY{a{i)HAxf" 


n=l 


C(Aa:) -f ln(a(t))-Fy^ 


2n (2n-|-l)! 

-lY{a{t)HAxY^ 

n=l 


ln(a(2")) ln(a(2)) 

3' 


2n (2n-|-l)! 


+6 


(89) 
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Now, we need to evaluate the four double integrals in Eq. flH^ to get the VEV. Results of 
the integrals are given in Appendix IA3I Adding up Eqs. flAldj) . flAlSjl . flA17D and flAlSjl 
and substituting the sum into Eq. fl8^ . we hnd 


{Q\ dt"ipl{t'',x) dt ifl{t,x')\fl) 


H3D-8 Y^(P-1) 

23D-2^3D/2 


3< —- ln^(a(f')) + ln^(a(f'))ln(a(t)) 


ln^(a(t')) — 21n^(a(t'))ln(a(t)) • 


+ln^(a(t'))ln^(a(t))—2C(Aa;^ 

—2C^(Aa:) [ln^(a(t')) — 31n^(a(t'))ln(a(t))] +4C^(Aa:)ln(a(t'))ln(a(t)) 
-l)^{HAx)^^ 


31n^(a(t'))i 2 




^ n‘^{2n+l)\ 

ln(a(f')), , , ,, 

^ ) ^A n(a(t))- 


ln(a(t)) — ln(a(t')) 


51n^(a(t0) ^ |C 


ln"(a(t)) 


9 




2C(Aa:)- —)ln(a(t')) 


ln(a(t)) 


Sn 


31n^(a(t')) 


2n 


^£{Ax,n)\n{a{t))—J'{Ax,n)\n{a{t'))+G{Ax,n) 


—S{Ax,n) ln(a(t)) + ln(a(t')) 


ln^(a(t))^ln^(a(t')) 


8n 


Sn 


■G(Ax,n) 


„ oo p-1 

- 

2 “4^ ^ pn (p — n) (277,+1)! [2 (p—u)+1 ]! 


-l)P{HAx)^P 


ln(a(t)) —ln(a(t')) 


ln(a(t'))+C(Ax)- — 


ln(a(f')) C(Aa;) 1 


p 


p p^ 


C(Aa:)- 


2 p 


ln(a(f)) + ln(a(f')) 


C(Aa:) 1 

A o 

p 


1 ~ {-l)i{HAxf'i a2'?(t')(ln(a(t)) -ln(a(t')) + ^) -ln(a(f)) -ln(a(f')) - J 

iZ^Z^Z^ — 


-q=3p=2n=i g(g-p)[2(g-p)+l]!n(p-77)(2n+l)![2(p-n)+l]! 

where we dehned new functions of comoving spatial separation Ax and index n 

_n 2 iK N 3C(Aa:) 13 , ._C‘^{Ax) 2C{Ax) 15 


(90) 


T{Ax,n) =C^{Ax) 


+ 


and G{Ax,n) = - 


+ 


(91) 


n 477^ ’ 77 77^ 877 ^ 

Recall that the function C{Ax) is expressed as a power series in Eq. fl42D . The square and 
the cube of C(Ax) that appear in Eq. fl90|) can also be expressed as series expansions 


OO p—1 






^ 477(p—r7)(2r7+l)![2(p— 77 )+!]! ’ 


00 (j—1 p—1 


c=(au=-ei:e 




(92) 

(93) 


3(g-p)[2(j-p)+l]!ij(p-ij)(2n+l)!l2(p-n)+l]! 

We have completed evaluating each of the three VEVs that yields an 0{X^) correction to 
tree-order two-point correlation function fl70|) of the IR truncated scalar full held. Substitu¬ 
tion of those VEVs given in Eqs. (jSS]), fl86|) and fl90|l into Eq. fl68|l yields the total two-loop 
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correction as 


llln(a(t')) 


^2 h ^ d -10 r3p_i) fl!^M))ln^(a(t)) 


60 


6 


12 

ln2(a(t')) 


ln‘^(a(t))+C(Aa;) 
ln^(a(t)) —C^(Ax) 


llli/(a(t)) ^ ln^(a(tO) ^ 


12 

21n^(a(t')) 


4 ' 3 

-21n^(a(t'))ln(a(t)) +C^(Aa;) 


3/A„^41n(a(t')) 


ln(a(t)) 


*-E * n(2u + ir ln(a(i))-ln(a(i')) 


li/(a(t')) ln'(a(t')) 


n 


- ln(a(t)) +--- 


X 


n n" 


In (a(t)) + 


+ 


24 L 


ln"^(a(t)) —ln"^(a(t')) +- ln^(a(t)) — ln^(a(t')) 


21 n(a(t')) 1 

21n3(a(f)) 21n2(a(f')) 

/2C(Aa;) 

9 

1 2 

n 

3n 

V n? 

2n?^ 


ln(a(t'))+P(C,n) 


In (a(t)) In (a(t')) .r/^ n , / / xx , / //xx xi 

H- . ^ ^ -- J\f{C,n) ln(a(t)) + ln(a(t)) -V{C,n) 


An? 


An? 




oo p—1 

- 

2^^^^pn(p—?7.)(2?7,+l)![2(p—n)+l]! 


a^^Y'Y ln(a(t)) —ln(a(t')) ln(a(t'))+C(Aa;) — 


2p. 


, ln(a(t')) 

CiAx) 

H 

1 

p 

1 

p 

p2j 


C(Ax)- 


2p 


ln(a(t)) + ln(a(t')) 


1 C{Ax) 

, H 

\ p 

p2j 


12 


EEE- 


^=3 p=2n=l 

where we defined 


(-l)9(HAa;)2‘'|a2'?(t') 

ln(a(t))-ln(a(t')) + 2 

1 

7 

7 

g(g-p)[2(g- 

-p)-|-l]!n(p —n)(2n-|-l)![2(p —n)-|-l]! j 


, (94) 


Af{C,n) = 


C2(Aa;) C(Ax) 1 


n 


n? n? 


and V{C,n) = 


C^{Ax) 2C{Ax) 13 




n-" 


4n"^ 


(95) 


To get the result for a fixed physical distance, one can choose the comoving separation 
Ax = K/Ha{t') so that the physical distance a{t')Ax remains as a constant fraction K of 
the Hubble length, as in Eq. (1471) . 

Adding up tree, one and two-loop correlators flTOD . (1771) and fl94|) yields the two-point 
correlation function of the IR truncated full field at one and two-loop order, 

(n|(p(t, T') |H) ~ (H|(^(t, X)p{t\ f') l^)tree + (f^|<^(^, f') l^)l-loop 

-h(H|(p(t,f)<^(f,f')l^)2-loop • (96) 
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(To save space we don’t write the result explicitly.) Computation of this correlation function 
was one of the goals of this paper. We use it, in the next section, to obtain the variance of 
the IR truncated scalar at one and two-loop order. 


C. Fluctuations in Spacetime 


The variation Aip of the IR truncated scalar ip is dehned as the difference of fields at two 
events 

Aip{t,t'; X, x') = ip{t, x)—(p{t\x') . (97) 

Because the VEV of held ip is zero at any event, the VEV of variation Aip also vanishes. 
However, the variance 

t'] f, x’) = (ff I [A(p- (f2| A(p|f2)]^|f2) = {n\ (A(p)^|f2) 

= (r2|(p^(t,T)|n) —(f2|(p(t,x)(p(t', a;')|r2) —(r2|(p(t', a;')<^(t,T)|f2)-l-(f2|(p^(T, x')|r2) , (98) 


is not necessarily zero. The VEVs in the middle of the second line in Eq. (|98|1 are just 
the two-point correlation function of the IR truncated scalar obtained in Eq. fl96l) . The 
remaining VEVs can be obtained from this function by taking appropriate limits: 


(r2|(p^(f, a;)|r2) = (r2|(p^(f, T')|n) 

H^d-& 


H 


D-2 


Tp-l) 


ln(a(t)) 


(L)-l)22^-2 37r^ r2(f) 


hr (a(t))- 


2D-I^DI2 
y2 ^3D-10 Y^{D-1) 


(D _ 1)2 23/5-3 5 7r3W2 pa (A) 


ln^(a(t)) , (99) 


and 

(r2|<p^(T, a;')|r2) = (r2|<p^(T, a;)|r2) = hm(r2|(p^(t, x)|r2) , (100) 

which can directly be read off from Eq. ([99]). Thus, using Eqs. fl96|) . fl99|) and fllOOl) in 
Eq. fl98|) . one obtains the variance a\^{t,t';x,x'). In the next two sections, we consider its 
equal time and equal space limits, respectively. 


1. Fluctuations in Space 

The spatial variation of the IR truncated scalar is the difference of helds at equal time 
events Aip(t,x,x') = ip(t,x) — ip{t,x'). It derives from the fact that inhationary particle 
production is a random process, like everything else in quantum mechanics. At some places 
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the field strength is positive, at others it is negative; at some places the magnitude of the 
field strength increases at others it decreases. Although the expectation value of the spatial 
variation A(p{t,x,x') vanishes, the spatial variance a\^{t,x,x'), which can be obtained by 
taking the equal time limit of Eq. (j9H]), does not. We employ Eqs. (j^ and (IMjl in Eq. (EH]), 
to obtain 


X, x') = (f2| [A(p(t, X, x') — {Q\A(p{t, X, \Q) = (f2| [(p(f, x) — (p{t, x')] |f2) 
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Tree order spatial variance in Eq. (llOip . 
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for a fixed comoving separation (increasing physical distance) asymptotes to 
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(We used Eq. (14T|) in Eq. fll02p to get Eq. fll03p .i Time dependence of tree order spatial 
variance implies that the magnitude of spatial variation increases with time. This fact is a 
crucial argument of some cosmologists against putting any trust in results for expectation 
values. They argue that an expectation value, which is measure of average effect, is mislead¬ 
ing since the actual effect that would be perceived by any local observer is either an increase 
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or a decrease in the field strength as fiuctuations happen. Reference [27|, on the other hand, 
argues that, while there is certainly spatial and temporal variation in the actual effects, one 
can, under certain circumstances, roughly trust what expectation values are telling us. See 
also Ref. 0 which mentions the problems some cosmologists have with using expectation 
values and gives conditions under which we can trust them. 

To infer our result fllOip for a fixed physical distance we take comoving spatial separation 
Ax = K/Ha{t) as in Eq. fl47j) . use Eq. (j39]), and get 
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_ / \ \ 


, sin(R) 

ln(R)H-^-ci 
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(_l)n^2n^-2n(^) 

2n (2n+l)! 


, (104) 


2D-2„D/2 pg) 1 

which agrees with the tree order result given in Ref. j27| when 0 = 4. Equation fllOdp implies 
that the tree order variance freezes in to a not especially large value. If the physical distance 
is taken as (half) the Hubble length, then the quantity in the curly brackets approaches 
(0.02) 0.08 during inflation. The one-loop variance in Eq. fllOip . 
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where a = a(t)HAx, for a fixed physical distance a(t)Ax = K/H, becomes 
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2 

—loop 


H2D-6 r2(D-l) fR:2 / 5 RTA, , , 

(0-1)220-%^ r2(f) |t 2 V’ 2’ “T") ln(a(t)) 
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which grows negatively. If the physical distance is taken as (half) the Hubble length, then 
the quantity in the curly brackets approaches (1.03) 4.07 after fifty e-foldings. This means 
that the spatial variation does not grow as large as tree-order variance fll04p implies. Hence 
the VEV describes the behavior of the field better when the loop corrections are included. 
Next, we study the temporal fiuctuations of the scalar field in our model, considering the 
difference of fields at equal space events. 
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2. Fluctuations in Time 


The temporal variation is the difference of fields at equal space events A(p{t,t',x) = 
(p{t, x) x) whose VEV vanishes. The temporal variance can be obtained by taking the 

equal space limit of Eq. (|9H]). Using Eqs. fl^ and ([99]), we hnd 


t', x) = (n| [A(p(t, t', x) - (QjAip(t, t', f) |n)]^|n) = (f2| [(p(t, x)-ip{t', f)]^|f2) 



y 2 ^3D-io y^{D-1) 81n^(a(t)) 171n®(a(t')) 
(D-l)2 23^7r30/2 r3(f) 5 ^ 30 


ln(a(T))ln^(a(t)) 4 



The one-loop correction tells us that the temporal variation does not grow as fast as the 
tree-order result implies. Thus, the VEV describes the behavior of field better, when the 
quantum corrections are included. In the next section, we study the IR truncated scalar in 
the same model applying the SET instead of the QFT. We will see that it yields the same 
two-point correlation function and hence the same statistical properties we infer from it. 


IV. STOCHASTIC FIELD THEORETICAL ANALYSIS 


Although the free field in Eq. fl21l) involves only the IR modes, it is still a quantum field 


in the sense that the commutator of the field and its first time derivative does not vanish: 


x), (po{t, x')]= ^'^^Q{Ha{t)—H27rn)[u{t,k)u*{t,k)—u*{t,k)u{t,k)]e^^'^^ ^ ^(108) 


n40 



See Appendix [B] for the derivation of Eq. (I109p . Recall that we replaced the mode function in 
Eq. (j^ with its constant leading IR limit (j27jl to obtain the IR truncated free field To(t, x) 
in Eq. fl28|l . Thus, Eq. (11081) implies that the commutator l(fQ(t,x),(fQ(t,x')] vanishes when 
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the replacement is made. That means that (po(t, a;) behaves like a classical variable (j)Q{t,x) 
in which the associated annihilation and creation operators are just random variables. Such 
a held is called stochastic. The full stochastic held (f){t, x) in our model obeys 

1 




( 110 ) 


hrst obtained by Starobinsky in Ref. 
up to the held strength renormalization factor {1 + 6Z)~^. 


18| . Note that Eq. flllOp has the same form as Eq. 


A. Free Theory 


As discussed earlier, we take the IR limit of the mode function and retain only the 
super-horizon modes —except the zero-mode, to obtain IR truncated stochastic free held 
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where k = 2nHn as in 
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r(f) { 2 k)^ 


( 111 ) 




. The annihilation and creation operators and At i 
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271 ] to be two complex conjugate stochastic random variables 


Arj = ^ {an+il^H) and ^ {an-iPfi) 


(112) 


The real random variables a^t and (3^ are governed by a Gaussian probability distribution 
function 


P (0)^1, (3fi) ■ 
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g 2 n g 2^71 
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with mean zero and standard deviation one. Note that the factors of \/^ in the denomina¬ 
tors exist for each mode n. 

Now, let us apply the stochastic formalism to compute the two-point correlation function 
of the IR truncated stochastic free held 
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where k' = 27iHrh. When m = n we have k' = k. Therefore, we obtain 
r2(D-l) 77^-2 




Using the fact that 
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where 9 is the angle between k = 2TTHn and Ax = x — x', and making the integral approxi¬ 
mation to the discrete mode snm over n as in Eqs. fl5^ - fl5^ we find 
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where we nsed Eqs. (|3^ and (HTll in the last eqnality. Let ns note here that the eqnal space 
and equal spacetime limits of Eq. (IllSh yield 

^ e(Ha(t)-k)ejHa{t')-k) , (ug) 


n 


n^O n^O 

where t' <t. Using Eqs. flll7p - flll8l) in Eq. (11151) . we find 


r(77-i) 


(l)o{t,x)(t)o{t\x 

770-2 r(D-i) 


2D-1t,D/2 

77^-2 r(D-l) 


./ //X X sin(a(f')hfAx) , sin(77Ax) 

ci{a{t')HAx) - \ \ [ -^-ci(77Ax)-' ^ ’ 


2D-i^d/ 2 rg) 


C(Ax)-|-ln(a(f')) 


a{t')HAx HAx 

'' [-lY{a{t')HAxf^ 

n=l 


( 120 ) 

( 121 ) 


2n (2?7,-|-1)! 

respectively. Results fll20p - fll2ip are the same as results fH51) - flT5p obtained in the context 
of QFT. Hence, 

(|0o(U x'jj = {kl\(po{t, x)(po{t\ f') |H) . (122) 

Because it will be needed in Sec. II V Cl let us compute the two-point correlation function 
for the hrst time derivative of the IR truncated stochastic free held. 
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which is 
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When rh = n we have k' = k, hence we obtain 
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Employing Eq. flllbp . making the integral approximation to the snm as in Eqs. flS4D - 
we find 
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This resnlt is the same as the one in Eq. flFT)) . Thns, snbstitnting Eq. fll25p into Eq. fll24p 
we obtain 
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which is the same resnlt obtained in Eq. fl58P applying QFT. The eqnal space limit of 
correlation fll26p is the limit given in Eq. fl59p . hence is a stochastic white noise. The 
limit enters the Fokker-Planck equation as a source term in the last section. In the next 
section, we study the self-interacting theory in the context of SET. 


B. Interacting Theory 


Here also we consider V{(j)) = ^0^ theory in D spacetime dimensions. Equation flllOp 


implies for the stochastic scalar 
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Iterating this equation twice for the IR truncated field 0 yields 
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Analogous to Eq. fl 68 |) . the two-point correlation function of the IR truncated stochastic full 
held is 


f')) = (0o(f, ^)0o(f', X') 
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where 0 < t'" <t" <t, 0<i<i<t' and t' < t. The hrst term is the tree-order contribution 
calculated in Eq. fll2ip . The 0{X) contribution involves two expectation values each of which 
has four free helds multiplied together. Recall that, in the quantum held theoretical analysis 
of Sec. im we reduced the computation to that of the topologically inequivalent pairings of 
free helds. These pairings originated directly from the commutator algebra [A^, At] = 
of the annihilation and the creation operators and the fact that = 0 which exclude any 

higher order groupings of the free helds. In the stochastic analysis, on the other hand, the 
Gaussian integrals over the random variables afi and [dfi do not vanish when an even power 
of the variables and jdfi are considered. Thus, one needs to exclude those combinations 
to allow only the pairings of the a^s and of the (das. This can be achieved by a projection 
operator V. Let us now compute the 0{X) contribution in Eq. (11291) . 


1. 0{X) Contribution 


The hrst of two expectation values which contribute at 0{X) in Eq. (11291) is 
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Inserting free held expansion (lllip . for each 0o, into Eq. (I130p . evaluating the Gaussian 
integrals using ^ and projecting out quartic combinations 
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of the random variables yields 




H 


2D-4 r4 


r^p-1). 


2iD-^Tr 


4~2D-2 


r‘(f) 


X 


dij2 

n^O 


Q{Ha{t) — k) Q{Ha{i) — k) 


n 


D-l 


cos 


k-{x—x'yj 


Q{Ha{i)-k’) 


m^O 


m 


D-l 


(131) 


where t<t' <t. At this point, we employ sums (11181) and flllOp in Eq. fll3ip and hnd 
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1=1 

This result is the same as the one obtained in Eq. fl72l) . Performing the remaining temporal 
integration yields the result given in Eq. (ITSP in the context of QFT. Thus, 


(j)o{t,x)l di(j)o{i,x')) = {kl\(po{t,x) I di(pl{t,x' 


The remaining expectation value which contributes at 0{X) in Eq. fll29p . 
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can be computed similarly. For each free held in Eq. (I134p we insert expansion fillip and 
evaluate the Gaussian integrals projecting out quartic combinations of the random variables, 
to obtain 
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We then substitute the sum computed in Eq flllOp into Eq. fll35p . 
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Because we have 0 <t"<t and t' < t, we decompose the integral into two parts: J^dt" = 
dt" + J^/dt" to evaluate the remaining sum in the integrand. In the hrst integral on the 
right t” <t', whereas in the second t' <t''. Employing Eq. 01161) in Eq. 0136p appropriately, 
we get 
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This result is the same as the one obtained in Eq. (175|) . We then perform the remaining 
temporal integrations and hnd the result given in Eq. fl76|l in the context of QFT. Thus, 


dfcplit”, x)(j)o{t', x')) = (fl| / dfiplit”, x)ipo{t', x')|f^) 


(138) 


Equations 0133p and 0138p imply the equivalence of the C>(A)-stochastic and the one-loop 
QFT correlators 


U{t,x)(j){t',x')\ = {Q\ip{t, x)ip{t', X ')\Q) i_ioop , (139) 

' / 0{A) 

where the latter was obtained in Eq. fITTj) . In the next section, we show that the same 
equivalence also holds for the C>(A^)-stochastic and the two-loop QFT correlators. 


2. C>(A^) Contribution 


Three expectation values contribute at 0{\^) in the two-point correlation function of the 
stochastic scalar full held in Eq. (11291) . Each will be computed in this section. Let us start 
evaluating the hrst one: 
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where the ordering of the time parameters is 0 <t <t <t' <t. We follow the same steps as 
in the evaluation of Eq. fll30p or Eq. (11341) . For each free held we insert expansion fillip 























into Eq. fll40p and evaluate the Gaussian integrals over the random variables projecting out 
quartic and hexic combinations. We hnd 
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where k” = 2nHr. The sums in Eq. (1141 p are evaluated in Eqs. (IllSp and (I119p . When their 
outcomes are substituted, the expectation value becomes 
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This result is the same as the one in Eq. fl82p . We, after integration, obtain the result given 
in Eq. ([83]). Thus, we conclude that 


(j)o{t,x) dt(j)Q{t,x') dt(j)'^{t,x')) = {kl\ipo{t,x) dt(pQ{t,x') dt(pQ{t,x')\kl) . (143) 


Next, we consider the second expectation value which contributes at 0{X‘^) in Eq. fll29p . 

































33 


Proceeding as before, we perform the Gaussian integrals projecting out quartic and hexic 
combinations of the random variables to obtain 
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Now, we use Eqs. (Illbh and (lllbh to rewrite Eq. (11451) as 

H2D-i r2(D-l) 


rt _ pt" 

+" J .2 (+n / j+m i 3 f+m 


1^ dt"0o(t", x)J^ dt”’(j)o{t”’, x)(j)o{t\ f'y - ^ 2 D- 2 ^D 


I dt" dt'" 

‘2 J ^ 0 0 


6(0o(t"x)0o(t',f'))lnl«(^"')) + 3(0o(t'"f)0o(t',x'))ln(a(t"))ln(a(r)) 


+6 0o(rx)0o(t',5^') ln^(«(^"')) 


(146) 


Recall that the stochastic and quantum formalisms are equivalent at tree-order, see Eq. (11221) . 
Thus, Eq. (11461) is identical to Eq. (1551) . 
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The hnal expectation value that contributes at 0{X‘^) in Eq. (I129p is 
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We carry on evaluating the Gaussian integrals projecting out quartic and hexic combinations 
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of the random variables and get, 
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At this stage we employ Eqs. flllSp and flllQp in Eq. fll49p . 
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The hnal step is to make use of the identity (I122p which implies that Eq. (IL^Op is identical 
to Eq. 


rt pf 

4-"Jki{+" I A+Ai(+ 


pt pf 
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(151) 


Equivalence of the individual stochastic expectation values to their corresponding quantum 
analogs, Eqs. fll43p . fll47p and fllSip . yields 

x)(j){t', x'jj = x)(:p{t', x')\n) 2 -loop ■ (152) 


Thus, Eqs. fll22p . fll39p and fll52p show that the stochastic held theory gives the same 
two-point correlation function fl96l) of the IR truncated scalar obtained applying QFT 


(^{t,x)(l)(t',x'fj = {fl\ip(t,x)(p(t',x')\Q) . (153) 

The variance of the stochastic held ^(t,x) dehned analogously to Eq. (j98p also yields the 
same quantum held theoretical result whose equal time and equal space limits are given in 
Eqs. fllOlIl and (11071) . respectively. In the next section, we give a stochastic check of the 
two-point correlation function of the IR truncated scalar in the equal spacetime limit. 
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C. A check of the via the Fokker-Planck equation 


Recall that we computed the two-point correlation function {VL\(p{t,x)(p{t\x')\VL) in 
Sec. mil using QFT in D-dimensions. We have just shown in Sec. IIVI that the same re¬ 
sult is obtained when SFT is used. In this section, we want to give another stochastic check 
of our result for the correlation function. 

The equal spacetime limits of Eqs. fl96|l and fll53p imply 


(r2|(p^(t, F)|r2) = x) 
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(154) 


5(F)-1)2 230-37r3W2 r3(f) 

One can also calculate the expectation values of functionals of a stochastic field 0 as 

/F[^{t,x)]\= [ duQ{t,u)F{u) , (155) 


where the probability density g{t, 0) is a solution of the Fokker-Planck equation 
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The expectation value we want to check in Eq. (11541) is of the form 


0^"'(f,x)^= f dujg{t,uj)uj‘^"' , 


(157) 


in F-dimensions with n = l. The computation of Eq. fll57|) was done in Ref. mj for D = A. 
Here, we need to do the computation for an arbitrary D. Taking the derivative of both sides 
of Eq. fll57p with respect to t, using Eqs. fll26p and fll56p . and then integrating the result 
by parts, we hnd 
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and choosing V{(p) = in Eq. (11581) we obtain 
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da 
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31| . So does its solntion 


This eqnation has the same form obtained in Ref. 

' ^Jd- 2 ^) \ = (2n-l)!!a4l-A-(n+l)a2 + P —(35n3+170n2+225n+74)a^ + C>(P) 
±1. j I - 2 280 


which yields, 




( 161 ) 


for n = 1. Employing Eq. fll59p in Eq. fll6ip reprodnces Eq. fll54p . The eqnal spacetime 
limit of onr resnlt fl96l) is indeed correct. 


V. CONCLUSIONS 

We considered infrared trnncated massless minimally conpled scalar field (p(t,x) with 
a qnartic self-interaction on a H-dimensional locally de Sitter backgronnd of an inflating 
spacetime. In Sec. IIIIl we compnted the two-point correlation fnnction {Q\(f{t,x)(f{t',x')\Q) 
of the scalar at one and two-loop order applying qnantnm field theory. The tree-order 
correlator at a fixed comoving separation, that is at increasing physical distance, approaches 
a constant nonzero valne. It grows linearly with comoving time at fixed physical distance. 
This is becanse more and more particles are created ont of vacunm during inflation which 
increases local field strength. The dominant (one-loop) quantum correction, however, grows 
negatively at a fixed comoving separation and at a fixed physical distance. 

We used our correlation function to compute the variance 

t'-, f, x') = (U| [A(^-(U| A(^|U)]^|U) = (U| (A(^)^|U) 

= a;)|r2) —(f2|(^(f, x)(f{t\ x')|r2) — (r2|(^(T, x')(f{t, x)|r2)-(-(r2|(^^(f', T')|r2) , 

and obtained corrections for the variance, in space and time, of one and two-loop correlators. 
Time dependence of the tree-order variance implies that the magnitude of spatial variation 
A(^(f, X, x') x)—(f{t, x') increases with time and freezes in, for a fixed physical distance, 

to a not especially large value. This spatial variation is the main reason for some cosmologists 
to argue that using expectation values gives a misleading description of a physical processes 
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because fluctuations about the mean value change the entire picture. A local observer 
detects either an increase or a decrease in the held strength as the huctuations take place. 
The expectation value measures merely the average effect. The authors of Ref. 27|, on 
the other hand, argue that one can roughly trust implications of expectation values under 
certain circumstances, despite the spatial and temporal variation in the real effects. The key 
issue in this paper is that whether the tree-order comparisons between expectation values 
and variances found in Ref. become better or worse when loop corrections are included. 
We found out that the one-loop variance grows negatively for hxed physical distance in our 
self-interacting model. This means that the spatial variance does not grow as large as the 
tree-order variance indicates. Thus, the contrast between the expectation value and variance 
decreases when we include the loop corrections. 

In Sec. HVl we examined the model by repeating the computations of Sec. [nil applying 
a stochastic field theory where the annihilation and creation operators of the free quantum 
held are considered to be two complex conjugate stochastic random variables. The real and 
imaginary parts of the random variables are governed by a Gaussian probability distribution 
function. We computed the two-point correlation function of the IR truncated stochastic 
held in our model. The result we found is the same as the one in Sec. m Thus, the 
analogously dehned variance in the context of stochastic held theory is in perfect agreement 
with the variance obtained using quantum held theory. Finally, we presented a check of our 
computation for the correlation function. Expectation values of functionals of a stochastic 
held can be calculated by integrating the functional weighted with a probability density 
function which satishes a Fokker-Planck equation. Employing this technique, we computed 
the expectation value of the stochastic held squared (0^(t,x)) in D-dimensions. The result 
is what our two-point correlation function yields in the equal spacetime limit. 
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Appendix A: Computing the two-loop contribution to the two-point correlation 

function {Q\(p{t,x)(p{t',x')\Q) 
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Two-loop contribution to the two-point correlation function of the IR truncated scalar 
full held consists of three terms each of which is proportional to a VEV; see Eq. (jMll . The 
VEVs are computed in Sec. IIIIB 31 In this Appendix, we give the details of the computation 
for each of the VEVs in sections lA 11 lA 21 and lA 31 respectively. 


1. Computing the VEV (n|(/3o(t, ®)/g dt <^o(b ® ^01^) 


In Eq. fl82l) computation this VEV is reduced to the evaluations of three double integrals. 
The hrst double integral is 
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The second double integral in Eq. flS^ is 
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The third double integral in Eq. 
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2. Computing the VEV {Vl\^^dt"ipQ{t", x)dt'"ipQ{t"\x)kpQ{t\x')\Q) 


In Eq. (18^ computation this VEV is reduced to the evaluations of three double integrals. 
Let us evaluate the hrst one. Because t' < t, we can break up the double integral into two 
(to be able to use Eq. (H5|) h 


pt pt" pt' pt" pt pt" 

/ dt" dt'"= / dt" dt'"+ dt" df . 

'o Jo Jo Jo Jt' Jo 


(A4) 


In the hrst double integral on the right we have 0<t'' <t', whereas t' <t" <t in the second 
double integral. Thus, using Eqs. (H5ll and (1A411 we can write the hrst double integral in 
Eq. fl85ll as 
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Let us now evaluate the second double integral in Eq. ([85]). To break up the double 
integral appropriately, consider the three cases that we may have depending on t': (i) 0 < 
t' <t'” <t” <t, (ii) 0 <t'" <t' <t” <t and (hi) 0<t'” <t” <t' <t. In cases (i) and (ii) we have 
t'<t", therefore, the second double integral in Eq. flA4p can be written as 
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In case (hi), on the other hand, we have t" <t'. Thus, in this case, the second double integral 
in Eq. flA4p can be written as 
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Note that the last equality in Eq. flASp is the same as the one in Eq. flA7p . Combining 
Eqs. flA4ll - flA8l) . we conclude that we can break up the double integral in Eq. flA4l) . in all 
possible cases, as 


ft pf 

Jt / Ajtt 


ff pf 

Jt / Ajtt 


ft pf 

Jt / Ajtt 


J pf 

Jt / Ajtt 


/ dt" dt'"= / dt" dt"'+ dt" dt'"+ dt" dt'". 

lo Jo Jo Jo Jf Jo Jf Jf 


(A9) 


























40 


Thus, using Eqs. (I45|) and flAOp the second double integral in Eq. becomes 
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The third double integral in Eq. fl55]l is obtained similarly. Using Eqs. fITS]) and flAOD we 
obtain 
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3. Computing the VEV {Q,\J^dt"YoY'Y)jQdiip'Q{i,x')\^) 


In Eq. (I89|) computation this VEV is reduced to the evaluations of four double integrals. 
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The first double integral in Eq. fl89|) yields, 
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The fourth double integral in Eq. flH^ yields 
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Appendix B: Computing the commutator [pQ{t,x),(po{t,x')] 


Here, we compute the commutator of the free scalar field given in Eq. (|2T|) and its hrst 
time derivative 
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where 


j dyy^ ^sm{y) = {D-3)y^ ^sm{y)-y^ ^ cos {y)-(p-3) (D-4)jdyy^ ^ sm{y) . (B4) 

Finally, using Eqs. flB2p . flB3p and flB4p in Eq. (IBip we get the commutator as 
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